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ABSTRACT 

We investigate the utility of the asymptotic giant branch (AGB) and the red giant branch (RGB) 
as probes of the star formation history (SFH) of the nearby [D = 2.5 Mpc) dwarf irregular galaxy, 
KKH 98. Near-infrared (IR) Keck Laser Guide Star Adaptive Optics (AO) images resolve 592 IR 
bright stars reaching over 1 magnitude below the Tip of the Red Giant Branch. Significantly deeper 
optical (F475W and F814W) Hubble Space Telescope images of the same field contain over 2500 stars, 
reaching to the Red Clump and the Main Sequence turn-off for 0.5 Gyr old populations. Compared 
to the optical color magnitude diagram (CMD), the near-IR CMD shows significantly tighter AGB 
sequences, providing a good probe of the intermediate age (0.5 - 5 Gyr) populations. We match 
observed CMDs with stellar evolution models to recover the SFH of KKH 98. On average, the galaxy 
has experienced relatively constant low-level star formation (5 x 10~ 4 M© yr -1 ) for much of cosmic 
time. Except for the youngest main sequence populations (age < 0.1 Gyr), which are typically fainter 
than the AO data flux limit, the SFH estimated from the the 592 IR bright stars is a reasonable match 
to that derived from the much larger optical data set. Differences between the optical and IR derived 
SFHs for 0.1-1 Gyr populations suggest that current stellar evolution models may be over-producing 
the AGB by as much as a factor of three in this galaxy. At the depth of the AO data, the IR luminous 
stars are not crowded. Therefore these techniques can potentially be used to determine the stellar 
populations of galaxies at significantly further distances. 

Subject headings: galaxies: stellar content — galaxies: irregular — stars: AGB and post- AGB — stars: 
color-magnitude diagrams (HR diagram) — instrumentation: adaptive optics 



1. INTRODUCTION 

Star formation histories (SFHs) are a key measure- 
ment for understanding galaxy evolution in a A Cold 
Dark Matter universe. Two main approaches have been 
used to constrain the SFHs of galaxies. One is to 
track the properties of different galaxy classes over cos- 
mic time, such as star formation rat e as a function of 
stellar mass (e.g. INoeske et al.l [20071 ) . The other ap- 
proach is to observe galaxies in the nearby universe 
and estim ate their SFHs based on their current prop- 
erties (e.g. iHeavens et ail 120041) . Results from both ap- 
proaches suggest a "downsizing" in star formation activ- 
ity whereby the most massive galaxies form the bulk of 
their stars early (before z ~ 1), while less massive, late 
type galaxies form the m aj ority of their stars l ater (e.g. 
iBrinchmann fc Ellisll2000T) . INoeske et all ([20071) suggests 
that both the onset and duration of star formation may 
be correlated with mass, so called "staged" star forma- 
tion. For gas-rich dwarf systems, this scenario should 
result in relatively constant low-level star formation over 
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much of a galaxy's lifetime. In this paper we derive the 
stellar populations for the local dwarf irregular galaxy 
KKH 98 as determined by optical/IR color- magnitude 
diagrams (CMDs) of the individual member stars. 

For the most nearby galaxies (d < 4 Mpc), color- 
magnitude diagrams of individual stars can be used 
to estimate SFHs. Optical Hubble Space Tele- 
scope (HST) imaging has been used extensively 
to study the stellar population s of ne arby galax- 
ies (e. g. ISaraiedini fc Jablonkal 120051: IWeisz et all 
2001 iGogarten et aU ^ I2009t IMcQuinn et all 120091: 



Williams et all I2009blfa1: iReikuba et all I2009D . For 



instance, deep HST and ground based imaging of Local 
Group dwarf galaxies has revealed a subset of dwarf 
spheroidals and ellipticals with extended star formation 
histories and star formation epi s odes as recent as 1 
Gyr ago (e.g. lAparicio et aTll200iT: IDolphin et all [20051) . 
The ACS Nearby Ga laxy Survey Treasury (ANGST; 
iDalcanton et all I2009D has expanded the Local Group 
efforts to include a large sample of more distant galaxies 
out to a distance, d < 4 Mpc. Using the HST Advanced 
Camera for Surveys (ACS) and archival HST imaging, 
ANGST has released photometry for deep two-color 
optical images of nearly 70 nearby galaxies. Matching 
the CMDs of these galaxies to model CMDs allows 
reconstruction of the star formation histories of a wide 
variety of galaxies in both group and field environments. 

Significantly less has been done to reconstruct the star 
formation histories from near-infrared (near-IR) stellar 
photometry. Several programs have used the SFHs of 
local group galaxies, as measured from optical data, to 
test stellar evolution models in the near-IR. For instance, 
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Fig. 1. — 40" X 40" images of KKH 98 in the F475W, F814W, and ftT'-band (top row; left, middle, and right respectively). Also shown 
are the three color composite (bottom left), and the full ACS field of the HST F814W-band image (bottom right). North is up and East to 
the left in all images. The spatial resolution of the Keck AO i^'-band image is a good match to the resolution of the optical HST images. 
The F475W and F814W band HST images reach depths of F475W~ 27.6 and F814W- 27.0 (12% photometric uncertainties). Photometry 
of these images recovered 2539 stars within the AO field. At this depth there is significant crowding in the HST images. In contrast, the 
Keck AO i^'-band image probes the brightest 592 near-IR sources, and shows no significant crowding. The very bright star at the top of 
the AO image was used as a tip-tilt guide star for the Keck AO system. 



IGullieuszik et al, (2008j) use the star formation history 
of Leo II dSph, obtained from optical CMDs, to esti- 
mate the expected production of AGB stars. They find 
that the models over-predict the observed numbers of 
AGB stars by a s much as a factor of six. Likewise, 
Gullicuszik et al. (2007) demonstrate how the near-IR 
colors of red giant branch (RGB) stars can be used 
to estimate the metallicity distribution within a galaxy. 
Other studies have used the near-IR luminosity functions 
of the AGB and RG B populations to estimate SFHs of 
local group galaxies (jOlsen et al.ll2006tlDavidgil2009l ) 

In spite of the limited analysis to date, the near- 
IR has several features that make it attractive for 
probing SFH. First, for most galaxies the bul k of the 
stella r luminosity is emitted in the near-IR (jSawickil 
2002). Second, intermediate aged asymptotic giant 
branch (A GB) populations are both e xtremely luminous 
in the IR (lAaronson fc Mould! [19851: 1 Fro gel et all 119901: 
IMaraston et aLlfeOOfjl ). and lie in tight sequences in near- 
IR color-magnitude space. Third, the effects of red- 
dening from dust are significantly reduced at near-IR 
wavelengths. Fourth, RGB age/metallicity degeneracies 
are reduced for optical-IR colors compared to optical 



or IR data alone (jGullieuszik et all I2007T) . Fifth, fu- 
ture missions such as the James Webb Space Telescope 
(JWST) and adaptive optics on the Thirty Meter Tele- 
scope (TMT) will be IR optimized. These missions will 
potentially resolve individual stars in galaxies out to the 
Virgo Cluster (18 Mpc) and beyond. In order to take full 
advantage of these missions a clear understanding of the 
IR properties of the CMD are required. 

To test the efficacy of near-IR photometry for recon- 
structing the SFHs of nearby galaxies, we have begun a 
campaign to image ANGST galaxies with Keck Adap- 
tive Optics (AO) and HST WFC3. In this paper, we 
present results from Keck AO imaging of nearby dwarf 
irregular galaxy KKH 98. This galaxy was chosen be- 
cause of its extended SFH and suitability of nearby AO 
guide stars. We obtained if'-band (2.12 /j,m) Keck laser 
guide star (LGS) AO imaging of KKH 98, as part of 
the Center for Adaptive Opt ics Treasury Survey (CATS 
Melbourne et all 120051 12008f ). The spatial resolution of 
the AO image (~ 0.1") matched the spatial resolution of 
the HST F475W and F814W images from ANGST. 592 
individual IR bright stars (K < 23.5) were resolved in 
the if -band, while 2539 stars were resolved in the opti- 
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cal (F475W < 27.6). We used the photometry of these 
stars to estimate the distance, metallicity, and star for- 
mation history of the galaxy. 

This effort is the first to use resolved IR photometry of 
individual stars to reconstruct the complete star forma- 
tion history (SFH) of a galaxy beyond the Local Group. 
We compare the IR results to the well constrained, opti- 
cally derived SFH of KKH 98, providing an excellent op- 
portunity to test these techniques at longer wavelengths. 
By refining SFH recovery methods in the IR, this work 
will help prepare for programs with JWST and TMT on 
more distant galaxies. 

Throughout, we report Vega magnitudes and assume 
a A cold-dark-matter cosmology: a flat universe with 
Ho = 70 km/s/Mpc, Vt m = 0.3, and Vt K = 0.7. 

2. HST AND KECK AO OBSERVATIONS OF DWARF 
IRREGULAR GALAXY, KKH 98 

This study uses high-spatial resolution imaging data 
from both HST and Keck AO to estimate the stellar pop- 
ulations within the dw arf Irr galaxy KKH 9 8. KK H 98 
was first identified by IKarachentsev et all ([2001) as a 
low surface brightness galaxy in the Palomar Observa- 
tory Sky Survey II plates. It is a gas-rich, isolated dwarf 
irregular galaxy with an H I li ne width of W50 = 2 6 km/s 
(|Karachentsev et al.l l200ll ). IKarachentsev et al.l (|2002ft 
find a distance of 2.45 Mpc based on the tip of the red 
giant branch. This distance gives an absolute magni- 
tude of Mb = —10.78. Its properties are fairly typi- 
cal of the local dwarf irre gular population identified by 
IKarachentsev et al.l (|200lD . 

2.1. HST Observations 

HST observations of KKH 98 were obtained in both 
the F475W and F81 4W filters of ACS as part of the 
ANGST survey (see iDalcanton et al.1 120091 for the de- 
tails). Three observations of the galaxy were made in 
each filter. The combined images have an effective ex- 
posure time of 2265 seconds and 2280 seconds, respec- 
tively, and are shown in Figure [TJ The ACS field of view 
(202" x 202") encompasses the bulk of the galaxy (See 
Figure [J). 

2.2. Keck AO Observations 

On the night of November 24, 2007, we obtained K'- 
band (2.12 ^m) Keck LGS AO imaging of KKH 98. We 
used the wide-field (40" x 40") camera of the NIRC2 in- 
strument to achieve the largest field of view. With a pixel 
scale of 0.04"/pix, the wide-field camera under-samples 
the AO PSF (0.06" resolution), limiting the effective res- 
olution to ~ 0.1". However, as is discussed below, the 
loss of resolution did not affect our ability to measure the 
photometry of the uncrowded IR-bright stars in KKH 98. 

Individual exposures were one minute, with a dither 
applied after two successive exposures. During the ob- 
servations, the laser was fixed to the center of the field. 
Thus, as the telescope dithered, the laser dithered as 
well. This method was preferred over fixing the laser to 
the sky because of significantly lower overhead. In ad - 
dition, as was demonstrated in iSteinbring et all ([2008), 
this method provides a more uniform AO correction over 
the field. An R = 14.2 tip/tilt star was located ~ 20" 
north of the galaxy center (Figure [lj . The AO system 
used this star to track and correct for image motion, 



while the laser spot was used for higher order wavefront 
corrections. 

We follo wed the image red u ction method outlined pre- 
viously in Melbourne et al.l (|2005f ). Sky and flatfield 
frames were created from the individual science frames 
after masking sources. Frames were then flatfielded and 
sky subtracted. We corrected frames for known NIRC2 
camera distortions. We aligned images by centroiding on 
sources, and combined them with a clipped mean algo- 
rithm. The final reduced image has an effective exposure 
time of 45 minutes. Figure [1] shows the resulting AO im- 
age of KKH 98. It covers the central 40" x 40" region of 
the galaxy. 

We observed the UKIRT photometric standard star 
FS3 immediately after the galaxy observations. How- 
ever, the night was not photometric. Therefore we set 
the zeropoint for the AO data based on two bright 
unsaturated stars in the actual scien ce frame which 
have if s-band magnitudes from 2MASS (Skrutski e et al.l 
2006). Ultimately, we compared our photometry to 
stellar isochrones pr oduced in the if-band system 
(jBessell fe Brettl Il988f l . However the transformation to 
the if-band system is negligible, typically < 0.02 mag 
for red stars (includ ing the transformation from Ks to 
K ()Carpen tcr 200 1 ]) and the color terms to transform 
from K' to Ks (jWainscoat fe CowidH99l V Therefore, 
for simplicity, we heretofore report all AO magnitudes in 
the if-band Vega system. 

The standard star observations of FS3 were obtained 
with tip-tilt AO correction to stabilize image motion, but 
without higher order LGS AO correction from the laser 
spot. This strategy produced a high signal-to-noise ratio 
(S/N) observation of the profile of the wings of the AO 
PSF (which are primarily set by seeing and tip-tilt). We 
used these standard star observations in the reconstruc- 
tion of the AO PSF as outlined in the following section. 

3. PHOTOMETRY 

iWilliams etHI (|2009bt ) and lDalcanton et al.l (|2009| ) de- 
scribe the photometry of the HST ANGST images. Stars 
in the ACS images were identified a nd measured with 
the photometry routine DOLPHOT (|Dolphinl [2000h us- 
ing PSF fitting of individual sources with model HST 
PSFs from TinyTirrQ. Aperture corrections were mea- 
sured for uncrowded stars and then applied to PSF fitting 
results. Photometry lists were culled to remove extended 
objects and low S/N stars. Stars with S/N > 6 in both 
filters were included in the final catalogue. The matched 
F475W and F814W-band catalogue contains 11324 stars. 
Twelve percent photometric uncertainties were measured 
at F475W= 27.6 and F814W= 27.0. Cutting the cata- 
logue at these flux levels, and limiting to the field of view 
of the AO image results in a final set of 2539 stars. 

The Keck AO image is significantly shallower than the 
HST frames (Figure[T]). However, the very deep HST im- 
ages already pinpoint the locations of stars in KKH 98. 
We take advantage of the HST positions when measur- 
ing stars in the AO image, performing photometry on the 
sources already identified in the HST frames. The AO 
PSF varies significantly across the field; thus we choose to 
use aperture photometry with aperture corrections cal- 
culated across a grid of sub-regions in the image. Discus- 

6 http:/ /www. stsci.edu/software/tinytim/ 
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Fig. 2. — Bottom: the azimuthally averaged stellar intensity as a 
function of radius (scaled to unity at a radius of 0.4") for the stan- 
dard star FS 3 (solid line), and 3 stars in the AO image (black, 
green, and red diamonds). FS3 was observed with tip-tilt cor- 
rection on but no higher order AO correction. The shape of this 
profile is therefore set by the seeing and tip-tilt correction. This 
figure shows that for radii larger than ~ 0.2", the halo of the AO 
corrected PSF has a similar profile as the tip-tilt corrected stan- 
dard star PSF. We therefore can use the high S/N PSF halo profile 
from our standard star (FS3) to track the profile of the our AO 
PSF out to large radii. The profiles of the 3 stars from the AO 
frame, which have higher order wavefront corrections, are signif- 
icantly more peaked in their core than the standard star which 
only had tip-tilt correction. Top: the fraction of light within a 
given radius for one of the AO corrected stars. A vertical dashed 
line denoting a radius of 0.16" is shown. This radius is the size 
of the aperture used in our photometry, and it contains more than 
40% of the stellar flux. 

sions of the AO photometry — including the PSF, aper- 
ture corrections, and aperture photometry — are given 
below. 

3.1. AO PSF 

To perform accurate photometry, we need a model of 
the AO PSF. Unfortunately, the AO correction degrades 
spatially with distance from the laser spot (anisoplan- 
tism) and the tip-tilt star (anisokinetisism), producing 
a spatially varying PSF. We model the PSF using two 
components, (1) an AO corrected core with a FWHM of 
roughly the diffraction limit of the telescope, and (2) a 
halo with a profile width set by the atmospheric seeing 
and AO tip-tilt correction. While the shape of the PSF 
halo is relatively stable across the image, the core varies 
significantly with position. 

To get an accurate representation of the profile of the 
PSF core, we identify luminous, isolated stars within 25 
equal-area sub-regions of the image. After subtracting 
neighbors, we register and coadd the stars in each sub- 
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Fig. 3. — The aperture correction for an 0.16" aperture as a func- 
tion of separation from the center of the AO image. The correction 
increases linearly with distance to the image center. We fit all data 
points with a line, giving a unique correction for each star in the 
image. This fit also improves the aperture correction for stars on 
the outskirts of the image where there are few stars to make a good 
estimate of the PSF. 

region, and trace the PSF profile of the coadded image. 

Unfortunately, the PSF halo is spread over hundreds of 
pixels, and is much lower S/N than the diffraction limited 
core. Therefore its profile cannot be accurately traced 
even for the stacked PSFs described above. Instead we 
used images of the very bright flux standard star, FS 3, 
observed with tip-tilt AO correction on, and high order 
AO correction off. With high order correction off, light 
that would have been corrected into the PSF core by the 
AO system, remained in the halo. The resulting PSF has 
a very high S/N profile out to large distance (~ 2") from 
the PSF center, and is a good match to the halo of the 
AO PSF (Figure HJ). 

To produce a final PSF for each of the 25 sub-regions, 
we splice together the halo and core images, scaling the 
halo image to m ap smoothly onto t he im age of the core. 
As was done in iMelbourne et al.l ([2008D . we select the 
splice point to be the radius at which the flux in the 
azimuthally averaged core profile is 9 sigma above the sky 
noise, as this is a good match to bright AO corrected stars 
where the halo is easily observed. Figure [2] shows that 
the PSF profile is relatively independent of the particular 
choice of splice point, for splice points larger than several 
times the diffraction limit (e.g. ~ 0.2"). This figure also 
shows that the halo PSF profile from the standard star is 
a good match to the few AO corrected stars in the science 
frame that are bright enough to adequately detect the 
halo. 

3.2. AO Aperture Corrections 

The stellar photometry of the AO image was performed 
with small apertures on the high S/N core of the AO 
PSF. We chose a 4 pixel radius aperture (0.16"). This 
aperture is 50% larger than the typical FWHM of the 
stellar PSF, and maximizes the S/N of the aperture pho- 
tometry. We convert the small aperture photometry to 
total flux with aperture corrections calculated from the 
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Fig. 4. — Photometry of 500,000 artificial stars embedded at ran- 
dom locations in the AO image. The top panel shows the difference 
of input vs. measured magnitude while the bottom panel shows the 
standard deviation of the difference as a function of input magni- 
tude. The thick dashed line represents the standard deviation cal- 
culated from all of the stars. The thin lines represent the standard 
deviations of groups of stars with consecutively larger separations 
from the image center (6, 12, 18, and 24 arcsec from the image 
center, as marked). Because the PSF is tighter near the center of 
the AO image, the S/N is higher for the photometry done there. 
The photometry shows low photometric uncertainty (better than 
15%) for stars brighter than K = 22. 

25 model PSFs, one for each sub-region of the image. To- 
tal flux for the model PSFs was measured with a curve 
of growth technique out to 2". 

The aperture correction varies from 0.6 mags at the 
center of the field to ~ 1 magnitude near the edge, where 
a larger fraction of the light remains uncorrected. Figure 
[3] shows the variation of the aperture correction across 
the field. A linear fit to the data gives aperture correction 
as a function of distance from the image center for an 
arbitrary stellar position. The fit provides a more finely 
sampled measure of the aperture correction across the 
field and improves the estimates in the outskirts of the 
image where the PSF is poorly measured due to a lack 
of stars. We use this functional form for the aperture 
corrections in the final photometry. 

3.3. AO Photometry 

The deep HST images of KKH 98 identify the locations 
of stars in the galaxy. To translate the HST positions 
onto the Keck AO image, we first selected a set of 157 
bright stars that are easily visible in both the HST and 



AO images. We then used the IRAF CCTRAN program 
to transform the HST coordinates into the AO image 
coordinates, using, a quadratic transformation. Since 
the geometric distortions in the AO image were removed 
during image processing, this transformation worked well 
across the entire frame. 

We measured the if-band magnitude for each star in 
the HST star list that was located on the AO frame. 
For each star, bright neighbors were identified with Sex- 
tractor. Scaled model PSFs were used to subtract the 
bright neighbors using a PSF fitting technique. We then 
measured aperture photometry on the star of interest, 
applying the aperture correction appropriate for the ra- 
dial separation of the star from the center of the image 
(Figure [3]). A correction for the local sky was made, us- 
ing a median sky value measured in an annulus of radius 
of 1.0 — 1.2 arcsec. When making the sky measurement, 
stars identified in the sky annulus by Sextractor segmen- 
tation maps and the HST star lists were masked out. 

Two types of uncertainty dominate the error budget 
in these measurements, (1) the photometric uncertainty, 
set primarily by the Poisson sky noise, and (2) the un- 
certainty introduced by the spatially varying PSF and 
aperture corrections. The photometric Poisson uncer- 
tainty can be measured directly from the images. For 
each star, we measure the standard deviation of pixels 
in the nearby sky to estimate this uncertainty. The un- 
certainty introduced by the aperture corrections can be 
measured from the scatter in Figure [3] While the scatter 
is small near the center of the image, it increases signifi- 
cantly towards the edges of the image. 

To estimate any systematics in the measurement meth- 
ods, we also performed a Monte Carlo simulation, em- 
bedding 500,000 artificial stars into the image at random 
locations with sub-pixel sampling. Artificial stars ranged 
in magnitude from K = 17 to 24, but were all assumed 
to be identified in the corresponding HST image. We 
measured the photometry of the artificial stars in the ex- 
act same manner as the real stars in the HST catalogue. 
Figure |4] shows the standard deviation of the photome- 
try as a function of input i^-band magnitude. We find 
to K = 22, the photometry is better than 0.2 mags. At 
K = 23.5 the photometric uncertainty is about a magni- 
tude. 

An additional 22 IR-bright stars were identified in the 
AO image that were not included in the final HST pho- 
tometry lists, usually because they fell below a blending 
or S/N threshold in the F475W image. We found that 
these stars do not alter the conclusions of this paper and 
therefore do not include them in the final analysis. 

4. RESULTS 

Figured shows the F475W - F814W vs. .F814W CMD 
and the F814W-ii~ vs. K CMD for KKH 98. Both 
CMDs show a tight red giant sequence with a tip at 
F814W= 23.2±0.1 and K ~ 21.2±0.1 (based on a visual 
fit to the CMD). The optical CMD exhibits a significant 
blue upper main sequence (F475W — F814W ~ 0.0) sug- 
gesting the presence of very young stars in this galaxy. 
The young red helium burning sequence is poorly popu- 
lated, although there do appear to be some young blue 
core helium burning stars. The main sequence is not ob- 
vious in the ¥%\AW-K CMD, because of the flux limits 
of the Jf-band data. 
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Fig. 5. — Apparent (uncorrected for reddening) color-magnitude diagrams of stars in KKH 98. F475W — F814W color is plotted on left, 
and F814W— K color is plotted on the right. Suspected AGB stars are shown as filled red circles. While both the optical and IR CMDs 
show a tight red giant branch, the IR color shows a much tighter AGB sequence, resulting in a powerful age indicator for intermediate aged 
populations. The 50% completeness limits (blue curve), and the cuts applied to the data, for modeling purposes (green lines) are shown 
for the optical data. The data cuts applied to the IR data are also shown. Typical photometric uncertainties at the faint end of each CMD 
are shown. At the bright end the uncertainties are roughly the size of the data-points. Data at the faint end of the IR CMD (K < 23) 
have photometric uncertainties greater than 0.6 magnitudes. Thus the large scatter in the IR CMD at faint magnitudes is not a feature of 
stellar evolution but rather photometric error. 



In the IR CMD, a well-defined asymptotic giant branch 
(AGB) emerges from the tip of the red giant branch 
(TRGB) and extends to K = 18.5. Compared with 
the optical CMD, the AGB stars in the F814W--K" 
CMD form a much tighter sequence. Tighter AGB se- 
quences are expected at redder wavebands for several 
reasons: (1) lower amounts of self-extinction from dusty 
AGB en velopes, (2) a natter co lor vs. stellar T e ff rela- 
tion (e.g. lGullieuszik et al.l l2007l. and (3) the thermally- 
pulsin g AGB stars have much smaller pulsation ampli- 
tudes |Whitelock et all 120061 : iMarigo et afl 120081 ) . Be- 
cause the AGB produces tight sequences in near-IR color- 
magnitude space, the AGB should prove a powerful in- 
dicator of star formation rate for intermediate aged pop- 
ulations, i.e. older than ~ 0.1 Gyr. 

At the faint end of the IR CMD (i.e. K > 23), pho- 
tometric uncertainties in the isT-band produce significant 
scatter in the photometry. Thus the large scatter to red 
colors in Figure [5] (and Figure [5]) is not a physical fea- 
ture in the data (such as the red clump), just photometric 
noise. 

Because of the spatially varying AO PSF, we investi- 
gated the possibility that the AO photometry might vary 
systematically across the image. However, CMDs gener- 
ated from multiple sub regions across the image show no 
systematic deviations. 

4.1. Modeling the Color-Magnitude Diagrams 

We estimate the star formation history (SFH) of 
KKH 98 by comparing the observed CMD s with CMDs 
produced from the stellar isochrones of iGirardi et al.l 



(2002), with updated AGB models from IMarigo et al.l 
(2008). These isochrones were used together with the 
bolometric correctio ns and color transformations from 
IGirardi et al.l (|2008|) . Sample isochrones, overlayed on 
the CMDs, are shown in Figure [Sj Metal poor (l/10th 
solar) isochrones (top row Figure O span the photomet- 
ric data. In contrast, more metal rich isochrones (bottom 
row Figure[6]) tend to be redder than the data, suggesting 
that the stars in KKH 98 are met al poor. 

We use the MATCH package (jDolphinl I2002T) to fit 
the CMD for a star formation history, metallicity his- 
tory, dust reddening, and distance modulus (DM), for 
an assumed Salpeter IMF and binary fraction of 0.35. 
MATCH uses a Poisson Maximum Likelihood statistic 
to determine the best fit model and the 1 sigma uncer- 
tainties on that mode l. For a comple te description of the 
MATCH routine seelDolphinl (I2002T). 

IWilliams et"afl ((20071) " showed that the IMF choice does 
not affect the relative star formation rates (as a function 
of time), but can change the overall normalization. Our 
primary goals are to compare the SFH measured from 
the optical data to the SFH measured from the IR data, 
and to determine if the SFH is "bursty" or continuous. 
As these are relative measurements, neither will be sig- 
nificantly affected by choice of IMF or binary fraction. 

Additional inputs to MATCH include: (1) a range of 
acceptable distances and reddening values; (2) time step 
sizes; and (3) ranges and binning sizes in color-magnitude 
space. In general we, attempted to reproduce the choices 
made in previous ANGST work (|Williams et al.l [20071 
2009b). Distance ranges were allowed within ±0.5 mags 
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Fig. 6. — Same as Figure \E\ only now with model isochrones overplotted (from Marieo ct al. 2008). Top: metal poor isochrones, l/10th 
solar, spanning ages from 10 7 - 10 10 years (colored lines) nicely span the data. Bottom: more metal rich isochrones (dashed lines) with 
an age of 10 9 5 years, the age by which roughly 80% of the stars had formed, are also shown and are significantly redder than the data. 
This suggests that the stars of KKH 98 are metal poor. Note: the very red colors measured for some faint stars in the IR CMD (below the 
horizontal line) can be explained by photometric uncertainty (greater than 1 mag.), rather than some usual stellar population. 



of the dista nce modulus in the liter ature of DM = 
26.95 ±0.11 l|Karachentsev et all [20021 ) . Reddening was 
allowed to vary from Ay = 0.1 — 0.7, a reasonable range 
for a dwarf galaxy with foreground redde ning Ai = 0.24 
as measured in the iSchlegel et al.1 (|1998f ) Galactic dust 
maps. We use a series of 12 logarithmically spaced time- 
steps which allow us to explore the star formation history 
in both the recent (Myrs) and distant (Gyrs) past. The 
minimum sizes for color-magnitude bins are set by the 
uncertainties of the stellar isochrones and are roughly 
0.05 mags in color and 0.1 mags in magnitude. How- 
ever, when setting bin sizes, additional considerations 



include allowing for reasonable number statistics in each 
bin while providing the resolution necessary to resolve 
features in the CMDs. For the optical data we used 
color bins of size 0.1 mags and magnitude bins of size 
0.15 mags. We doubled the sizes of these bins in the 
NIR to offset the smaller number statistics in the NIR 
CMDs. Small changes in the bin sizes were not found to 
alter the results. 

We used the latest version of MATCH (version 2.3) 
which handles the isochrone color transformations for 
Carbon and Oxygen AGB stars independently. In older 
versions of MATCH, C/O ratios were not tracked and the 
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Fig. 7. — TOP: the optical CMD now displayed as a density 
Hess Diagram. The true data (top left), and the best-fit MATCH- 
derived model (top right) are shown. Also shown are the residual 
differences between the two (bottom left), and the residual differ- 
ence normalized by Poisson statistics in each bin (bottom right). 
In the residual maps, darker regions denote an under-production of 
model stars. Lighter bins denote an over-production of model stars. 
The residuals (bottom left) vary from +10 stars per bin (white) to 
-11 stars per bin (black). The normalized residuals (bottom right) 
vary from +8 (white) to -5 (black) and can be thought of as sigma 
deviations per bin. The models reproduce the data, except at the 
faint end, and for the AGB, where the models over-produce AGB 
stars by roughly a factor of three (F475W ~ 25.7 and F475W - 
F814W ~ 3, see Figure [9] for more details). 



color transformations produced AGB stars significantly 
redder (2-3 mags) than expected. The colour transfor- 
mations of C stars in the updated v ersion of MATCH 
were derived from iLoidl et ail ()2001[ ) . Because the in- 
formation content of the IR CMD is biased towards the 
AGB populations, it was important for us to use the most 
realistic prescriptions for AGB photometry. In the op- 
tical CMDs these considerations were less vital because 
the information content at young and intermediate ages 
primarily arises from the Main Sequence turn-offs. 

For the optical data, we chose to model the CMD with 
photometry brighter than the 12% uncertainty limits and 
50% completeness limits (F475W< 27.6 and F814W< 
27.0) to avoid fitting areas of the CMD with poor pho- 
tometry or large completeness corrections. These limits 
result in photometry for 2539 stars in the area overlap- 
ping the AO data. For the IR CMD, we chose stars with 
K < 23.5, the point at which the -ftf-band photometric 
uncertainty is about 1 magnitude. These cuts provide 
592 stars. While the photometric precision at the fainter 
end of K is poor, the long wavelength baseline between 
the F814W and fsT-bands means that there is still sta- 
tistically interesting information at these fainter limits. 
Table []] gives a summary of the MATCH derived results 
which are discussed in detail below. 

Figures [7] and [8] show Hess diagrams of the CMDs in 
the optical and IR (top left). Also shown are the best fit 
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Fig. 8. — Same as Figure [7] only now for the IR data. The 
residuals (bottom left) vary from +10 (white) stars to -9 stars 
(black), while the normalized residuals vary from +5 sigma to -5 
sigma. The fit to the AGB is better for the IR CMDs than the 
optical CMDs, although the models st ill o ver-predict the the AGB 
by roughly a factor of two (see Figure ITol for details). 

model CMDs derived by MATCH (top right), the resid- 
ual difference between the model and the data (bottom 
left), and the residuals scaled by the statistical signif- 
icance of bin (bottom right). The residual differences 
between the data and MATCH derived models are small 
for both the optical and IR CMDs, except for AGB stars, 
especially in the optical data where the models predict 
130 AGB stars compared to 42 observed (with the caveat 
that the AGB must be brighter than the TRGB, and 
redder than FA75W - F8UW = 1.75). In the IR, the 
models predict 92 AGB stars (brighter than the TRGB). 
We return to these discrepancies in the following section 
and elaborate in the discussion section. In contrast, the 
numbers of stars in other evolutionary stages, such as 
the red giant branch and main sequence are better re- 
produced by the models. Table [2] gives the measured 
and predicted numbers of stars in different regions of the 
KKH98CMD. 

4.1.1. Star Formation Histories 

The optical and IR derived SFHs of KKH 98 are plot- 
ted in Figure [TTJ For the bulk of cosmic time, the SFRs 
from the two CMDs agree to within the uncertinties. For 
ages older than 3 Gyrs, both data sets show evidence 
for a modest, roughly constant star formation rate of 
5 x 10 _4 A/ o yr-\ or 2 x 10" 4 A/ Q yr" 1 kpc~ 2 . From 
1 — 3 Gyrs ago, both show a drop in star formation to 
roughly negligible rates. More recently, both show re- 
newed star formation for ages < 0.5 Gyrs ago. However, 
in these youngest time bins, there are differences in the 
details. At the youngest ages (age < 0.1 Gyr), the op- 
tical photometry suggests ongoing star formation, while 
the IR data set does not. This difference is not alarming, 
as the near-IR data are not very sensitive to hot (blue) 
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Fig. 9. — A realization of the best fit optical model CMD shown 
broken into different time bins (blue diamonds). The observa- 
tional data are also shown (black points), with the AGB stars 
circled in red. The models shown that the main sequence is 
populated by young stars (log(age) < 8.8), whereas the red gi- 
ant branch is primarily populated by old stars (log(age) > 9.5). 
The model predicted numbers of AGB stars are given in the up- 
per right of each CMD. For AGB stars above the TRGB (and 
F475W - F814W > 1.5), the data show 42 stars, while the models 
predict 130 stars. The lower left plot shows the magnitude limits 
provided to MATCH (green lines). 

main sequence stars produced by recent star formation. 

Another difference is that the optical data set suggests 
a burst of star formation, at age ~ 0.5 Gyrs (Figure 
[TTj) . In contrast, the IR data suggest SFRs roughly a 
factor of four smaller at that time. The evidence for a 
burst in the optical data comes from stars at the well- 
populated main sequence turn off (MSTO) at the faint 
end of the optical CMD. This is shown in Figure |H] where 
a realization of the predicted model CMD divided into 
different time bins. At log(age) = 8.6 — 8.8 [yr], large 
numbers of stars are predicted at the MSTO. In contrast 
the IR data do not reach deep into the main sequence, 
and therefore the IR-derived SFH is not sensitive to the 
main sequence turn-off feature. However, the IR data 
should contain information on populations of that age 
from the AGB sequence (see Figure [T0|) . Therefore, if 
the burst is real, evidence for it should exist in the AGB 
population, assuming the AGB models are accurate. In 
the discussion section we examine possible explanations 
for the measured differences in the SFR at these ages. 
Here we discuss the effect of time binning. 

Time binning may play a role in the differences between 
the optical and IR derived SFHs. Logarithmic time bins 
have been selected to fully explore the SFH on both re- 
cent and older time scales. However, especially for the 
IR data, which has only 592 stars, the time binning may 
be more fine than the data warrant. The third panel of 
Figure [11] shows the IR results now binned with larger 
time steps. These bins were selected by first binning the 
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Fig. 10. — Same as Figurc[9] only now a realization of the best fit 
IR model CMD (blue diamonds) shown broken into different time 
bins. The model predicted numbers of AGB stars are given in the 
upper right of each CMD. For AGB stars above the TRGB, the 
data show 42 stars, while the models predict 93 stars. In this case 
we see that the bulk of the predicted AGB stars are old (log(age) 
< 9.5). The lower left plot shows the magnitude limits provided to 
MATCH (green lines). 

data on small time scales and then tracking the goodness 
of fit parameters as bins were dropped from the sample. 
If the fits did not change by more than the measured 
uncertainties, a dropped bin would be incorporated with 
its neighboring bin. Ultimately, this resulted in 4 bins. 
We caution, however, that this method is likely to miss 
the importance of bins with little star formation such as 
those from 1 — 3 Gyrs ago. For more on optimal binning 
techniques, see Williams et al. (2009 in press). The opti- 
cal data have been binned to match the IR bins. Binned 
this way, there are still some differences between the two 
measured star formation histories at the youngest ages. 

The SFR uncertainties shown in Figure [11] are a com- 
bination of (1) the systematic uncertainty between the 
models and the KKH 98 photometry, and (2) the ran- 
dom uncertainty associated with small number statistics 
of stars in any given color-magnitude bin. The first un- 
certainty is estimated by the MATCH routine from the 
maximum likelihood statistic as it attempts to fit models 
parameters within the adopted ranges. The second un- 
certainty is estimated by a Monte Carlo bootstrap sim- 
ulation of 100 realizations, where we rerun MATCH but 
alter the input photometry list, so that it is composed of 
random draws from the actual photometry list. The sys- 
tematic uncertainties are roughly equivalent for both the 
optical and IR data, and are similar in scale to the ran- 
dom uncertainties associated by small number statistics. 
The error bars shown in Figure [TT] are the two uncertain- 
ties added in quadrature. 

Surprisingly, the optical and IR model uncertainties 
are similar in scale despite the fact that the optical data 
contain contain a factor of 5 more stars. Several factors 
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Fig. 11.— The MATCH derived star formation history of KKH 98 
from the optical CMD (blue) and the IR CMD (red). At early 
times, both show a relatively constant star formation rate of 
5 X 1O" 4 M yr-\ or 2 X 10~ 4 M© yr" 1 kpc~ 2 . From 1-3 Gyrs 
ago, both show a significant drop in the star formation to negligible 
rates. More recently, (age < 1 Gyr) the SFR has roughly returned 
to its past average. However, the optical data suggest a possible 
burst of star formation at 0.5 Gyrs ago, not reproduced in the IR 
data. The IR CMD is also insensitive to the most recent star for- 
mation (i.e. age < 60 Myr). Uncertainty estimates are shown for 
each time bin, and combine in quadrature (1) the systematic un- 
certainty between the models and the data, and (2) the uncertainty 
associated with small number statistics in a given time bin. The 
bottom panel shows the data rebinned to longer time periods. 

conspire to lead to this result. First, the color- magnitude 
bins used in the IR analysis are twice as large as the 
optical analysis, and thus typically contain as many or 
more stars as the optical bins. Second, a lack of stars 
in a feature can be just as important to producing a 
significant result as the presence of stars. For example 
the lack of AGB stars in the data forces the IR models 
to select against a burst of star formation at 0.5 Gyrs 
ago (see Figure ITUl) , resulting in small numbers of AGB 
stars in both the models and data, and a low uncertainty 
assigned to this time bin. Whereas in the optical data, 
the main-sequence turn off suggests a burst at that time, 
and the lack of AGB stars acts to increase the uncertainty 
of the high SFR measured in this time bin. 

4.1.2. Metallicity 

Figure [T^] shows the MATCH measured metallicity of 
new stars at every epoch. Again results from both the 
optical (blue) and IR (red) CMDs are shown, however, 
in time bins where no star formation was measured, no 
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Fig. 12. — The MATCH derived metallicity of new stars at each 
epoch from the optical (blue asterix) and IR (red diamonds) CMDs. 
Both show a preference for metal poor stars from roughly l/30th 
to l/10th solar. There may be a slight increase in metallicity over 
time. For timebins with negligible star formation estimated in the 
IR CMDs, no metallicity is plotted. 

metallicities are shown (several IR time bins). Both sug- 
gest that the chemical history of KKH 98 has been dom- 
inated by low metallicity star formation, producing stars 
with typical [M/H] ~ —1.5 to —1.0. There may be a 
trend to more metal rich stars at more recent times. 
MATCH can also be run in a mode where the metal- 
licity is constrained to continually increase, as would be 
expected in a closed box galaxy model. When run in this 
mode, the MATCH results do not change appreciably. 

4.1.3. Reddening and Distance 

MATCH provides an estimate of the foreground dust 
obscuration and distance to KKH 98. By varying these 
two parameters, MATCH effectively normalizes the mod- 
els in color magnitude space so that they best match the 
observational data. Because different stellar evolution 
models will have different normalizations, the best fit 
reddening and distances reported by MATCH are model- 
dependent and therefore should not be considered mea- 
surements of the actual distance and foreground redden- 
ing. However, a comparison of the reddening and dis- 
tance results from the two CMDs presented in this paper 
will help inform whether consistent results are possible 
from the optical and IR data sets. 

For the optical data, MATCH infers a gl obal Ay — 
0.52 ± 0.06. Using the iCardelli et all ([l98l reddening 
law, this translates to an Aj = 0.23. To within the 
uncertainties, the IR data gives a similar dust absorp- 
tion of A v = 0.42 ± 0.08. Using the optical CMD, 
MATCH finds a distance modulus of DM = 27.07X0.09. 
The IR CMD produces a similar distance modulus of 
DM = 27.10 ± 0.07. To within the uncertainties, both 
estimates are similar to distan ces quoted in the literature 
- e.g. DM = 26 .95 X O.lKlKarachentsev et al.l l2002l) 
and DM = 27.023 (jDalcanton et all 120091 ) . More on the 
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Fig. 13. — The magnitude of the tip of the red giant branch 
(TRGB) as a function of metallicity and age in the F814W (top) 
and if -bands (bottom) as derived by the isochrones of Mari go et al.l 
l(200Sl . colored lines), the same isochrones used by MATCH. Also 
shown is the metallicity dependence on the i-C-band TRGB as esti- 
mated from old G alactic globular clusters (dashed black line from 
IValenti et aT][200i) . The MATCH-derived metallicity of KKH98 
of [M/H] = —1.4 is indicated by a hatched vertical region. For 
metal poor populations, M^g^Sr is roughly constant with metal- 
licity, and for the metallicity of KKH 98, AfJ^J^ = ~ 4 -0- For 
metal rich systems, the M^ RGB is roughly constant with metallic- 
ity. At the metallicity of KKH 98, the lValenti et al.l 1120041) relation 
suggests M^ RGB = —6.05, whereas the Marigo ct al. (2008) model 
isochrones suggest a K-b&nd tip about 0.2 mags fainter. The latter 
estimate gives a more consistent distance to KKH 98. 

distance to KKH 98 is given below. 

4.2. Estimating Distance from the Tip of the Red Giant 

Branch 

While the MATCH-derived distances arise from nor- 
malizing the entire CMD to models, specific features in 
the CMD can also be used to estimate the distance. One 
of the best distance indicators available in the KKH 98 
data set is the F814W-band magnitude of TRGB. Fig- 
ure [13] shows the F814W-band magnitude of the TRGB 
as a function of metalli city and age, using isochrones 
from IMarigo et al.l ( 20080 For old, metal-poor stellar 
populations the TRGB has an absolute magnitude of 
M Fsdw = -4.0 (Figure [T3]) . We estimate the apparent 
magnitude of the TRGB in KKH 98 to be at F814W= 
23.2±0.1. Correcting for reddening using Aj — 0.23, this 
gives a distance modulus (DM) of DM = 26.97 ± 0.10, 

7 Models from CMD 2.1 |http: / /stev.oapd. inaf.it /cmd] 



in good agreement with jKa rachcnts ev et"aLl ()2002l ) who 
found DM = 26.95 ± 0.11 from the J-band TRGB. 

The apparent Jf-band magnitude of the TRGB (K — 
21.2 ± 0.1 for KKH 98) can also be used to estimate 
the distance of KKH 98. In Galactic globular clusters 
with old stars, the if-band magnitude of t he TRGB has 
been s hown to be a function of metallicity (Vale ntTet al.l 
120041) . iGullieuszik et all (|2007t ) parameterize this empir- 
ical function as: 



m trgb = _ g 92 



0.62 [M/H], 



(1) 



Assuming the MATCH-estimated KKH 98 metallicity for 
old stars of [M/H] = -1.4, equation [1] gives M^ R ° B = 
—6.05. Given the apparent TRGB magnitude of K — 
21.2 ±0.1 and negligible dust absorption in the if-band, 
we calculate a KKH 98 distance modulus of DM — 
27.25 ± 0.10. This distance measurement is somewhat 
farther than the distance estimated from the apparent 
F814W-band magnitude of the TRGB. 

However, another estimate of the abso l ute if -band 
TRGB can be made from the IMarigo et all (|2008D mod- 
els shown in Figure[T3"l which are also the models used by 
MATCH. These models suggest that, at low metallicities 
and old ages, the if-band magnitude of the TRGB is 0.2 
magnit udes fainter than the value from Equation [TJ Us- 
ing the Mari gcTet al.l ([2008) models, we find a distance 
modulus of DM = 27.05 ± 0.10, which is in very good 
agreement with the DM derived from optical CMD. 

5. DISCUSSION 

Using high spatial resolution optical (HST) and near- 
IR (Keck AO) images, we have constructed CMDs of 
the stellar populations within dwarf irregular galaxy, 
KKH 98. From these data, we reconstructed the full 
star formation history of the galaxy, the first time this 
has been done in the IR for a galaxy beyond the Local 
Group. Below we compare the optical and IR results 
and demonstrate how these techniques may be useful for 
galaxies at larger distances. 

5.1. Comparison of Optical and IR Results 

Although the number of stars available for analysis in 
the IR was an order of magnitude smaller than for the 
deep HST optical data, the bulk of measured properties 
— including the broad trends in star formation history 
and metal enrichment — agree between the two data sets. 
There are some differences however. The most obvious is 
that the IR data set is not sensitive to the youngest stars 
(i.e. age < 0.1 Gyr), because the bluest main sequence 
stars emit the bulk of their radiation in the UV and op- 
tical. Therefore, HST optical observations are required 
to constrain the star formation histories at the youngest 
ages. 

There are also differences in the details of the star for- 
mation histories at ages of ~ 0.1 — 1 Gyrs ago, precisely 
when AGB stars are expected to dominate the informa- 
tion content of the IR CMD. For instance, 0.5 Gyrs ago 
the optical data suggest a burst with a star formation 
rate four times the prediction of the IR data. The evi- 
dence for this burst comes from the well-populated MS 
turn-off of the optical CMD at that age (Figure O- 

A burst of this age should produce AGB stars ob- 
servable in our optical and IR CMDs. Unfortunately, 
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Fig. 14. — Same as Figure 1101 only now the models (blue di- 
amonds) are created using the optically-derived SFH. The over- 
production of the AGB in the old time bins is basically identical 
to that for the IR derived star formation history. However, in 
this case, significant numbers of AGB stars are also predicted at 
younger ages. The young model AGB stars are not only over- 
produced, but especially for log(age) = 8.6 to 8.8, their colors are 
also typically 0.5 to 1 magnitude redder than the actual AGB pop- 
ulations in KKH 98. 

MATCH predicts that this burst should produce roughly 
the entire observed AGB population, which would only 
make sense if there were no AGB production at older 
times (Figures O • The over-production of AGB stars by 
the models is even more striking if we model the burst 
in the IR. Figure [14] shows a realization of the model IR 
CMD created with the optically-derived star formation 
history. Not only does the model predict large num- 
bers of AGB stars at 0.5 Gyrs ago, it also suggests that 
their colors should be 0.5 — 1 magnitude redder than the 
AGB populations seen in the data. As the SFH derived 
from the IR data relies heavily on the AGB sequences, 
it is now clear why the IR data under predicts the star 
formation in this time bin; there are few AGB stars in 
the KKH 98 data that match the luminous, red AGB 
stars predicted by the models for a burst of that age. 
Figure [H] also suggests that age/metallicity degeneracies 
are not producing the discrepancy; the predicted AGB 
stars are well separated from the observed AGB stars in 
color space. 

Not only do the models over-produce young metal poor 
AGB populations (and predict redder colors than are 
seen), they also over-produce old metal poor AGB pop- 
ulations (Figures [91 and IT0|) . For populations older the 
3 Gyrs, the models predict ~ 57 ± 8 (optical, assuming 
Poisson uncertainties) and ~ 79 ± 9 (IR) AGB stars in 
KKH 98. These predictions are both larger than the to- 
tal number of 42 observed AGB stars, which should pre- 
sumably include some younger AGB stars as well. Taken 
together, these results suggest issues with models of AGB 
stars especially for metal poor populations such as those 



found in KKH 98. 

The modelling of AGB depends crucially on the pre- 
scriptions for two difficult-to-model phenomena: mass 
loss and the efficiency of third dredge up episodes. To 
better address the uncertainties related to these pro- 
cesses, modellers resort to synthetic codes in which a 
few parameters are tuned to fit observ ational data, typi- 
cally from the Magellanic Clouds (e.g. iMarigo fc Girardl 
[200% . 

While synthetic AGB codes have successfully re- 
produced the numbers, luminositie s, and colors of 
AGB stars in the Magellanic Clouds (fcToni et al.l [19991; 
iNikolaev fc Weinberd[2000t IMarigo fc Girardill2007l ). the 
situation in other galaxies h a s been less encouraging. For 
instance, Gul lieuszik et al.l (|2008l ) showed that models 
calibrated on the LMC and SMC over-predict the AGB 
lifetimes of the old, metal poor dwarf spheroidal, Leo II, 
by as much as a factor of six. Girardi et al. (in prepa- 
ration) reaches a similar conclusion based on a larger 
sample of ANGST galaxies. 

Perhaps more relevant to the pre sent paper are the con- 
clusions from a similar work from lHeld et al.l (|2010D : In 
Leo I dSph, a metal poor galaxy containing young stel- 
lar populations (therefore more similar to KKH 98), they 
find that the carbon stars predicted by the models are 
about twice those observed. We may be seeing a similar 
issue in KKH 98, where we would need to roughly triple 
the AGB population in the data to match the models, 
and the difficulties with the AGB seem to be particularly 
problematic at younger ages. We will potentially be able 
to better constrain AGB models in this vital time range 
with observations of the AGB populations of a wide va- 
riety of galaxies, as will be possible with the ANGST 
WFC3 survey. 

5.2. Application of These Techniques in Other Systems 

In KKH 98, the IR bright stars are well separated, with 
typical separations of ~ 1". Assuming these separations 
are common within other galaxies, crowding should not 
limit the technique until much larger distances (up to 
10 Mpc for Keck). The primary limiting factor for fu- 
ture ground based observations with AO will therefore 
be sensitivity above the thermal background. Our study 
was completed with a 45 minute exposure. At larger dis- 
tances, significantly more exposure time, or better AO 
correction will be required to reach the depth of our 
study (~ 2 magnitudes below the TRGB). For instance, 
to reach comparable S/N in a galaxy twice as distant 
would require 16 times the total exposure time, or 12 
hours. There are, however, several ways in which this 
time requirement could be reduced. 

1) Improved AO performance. The observing condi- 
tions during this run were sub-optimal with only ~ 20% 
Strehl ratio (ratio of peak luminosity in the PSF to the- 
oretical maximum peak luminosity for a diffraction lim- 
ited image). With the Keck AO upgrade, typical Strehl 
ratio's of 30% or more are common. A higher Strehl ra- 
tio will dramatically increase the sensitivity above back- 
ground. In addition, if the Strehl ratio were more uni- 
form across the image larger numbers of faint stars would 
be well photometered towards the edges of the image. 
More uniform AO performance will be possible with up- 
coming Multi-Conjugate AO systems such as the one be- 
ing developed for Gemini south, which will have an iso- 
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planatic patch 4 times the size of the Keck AO system. 

2) At shorter wavelen gths, the thermal I R bac kground 
is reduced. In g-band. iMelbourne et al.l (|2007l) showed 
that it was possible to photometer an H = 24.0 [Vega] 
point source to a photometric precision of 17% in an hour 
long Keck AO exposure. The ff-band sensitivity is sig- 
nificantly deeper than the i^-band limit. The drawbacks 
for using this strategy are, (1) the AO system typically 
achieves smaller Strehl ratios at shorter wavelengths, and 
(2) there is a smaller wavelength lever-arm in the CMD 
when using ii-band (1.6 /jm) data compared with K- 
band (2.2 /im). Despite these limitations iJ-band may 
be a good strategy for galaxies at larger distances than 
KKH 98. 

3) With near-PR WFC3 observations from HST, both 
the thermal background and field of view issues can be 
solved. We expect the planned WFC3 observations of 
the ANGST sample to reach significantly deeper than 
our Keck AO observations. These observations will be 
made in both the H and J-bands. The roughly factor of 3 
loss in resolution should be relatively unimportant at the 
distances of the ANGST sample, where we have shown 
that the IR luminous stars are typically not crowding 
limited. In more distant galaxies, crowding will be a 
problem for HSTdX these wavelengths. 

4) Future instruments such as AO on a Thirty Meter 
Telescope or JWST should also overcome the limitation 
of our current data set. These telescopes will potentially 
resolve individual stars in galaxies out to the Virgo Clus- 
ter (18 Mpc) and beyond. 

6. CONCLUSIONS 

We obtained high spatial resolution if-band imaging 
of dwarf irregular galaxy KKH 98 using the Keck LGS 
AO system. The resolution of the AO images (~ 0.1") 
was a good match to the existing HST optical images. 
In addition, the AO PSFs had a significantly higher con- 
trast above the background compared with seeing limited 
images, resulting in 15% photometry at K ~ 22 [Vega] 
for our 45 minute exposure (si milar to the Keck A O pho- 
tometric uncertainty found bv iVacca et al.l I2007L in the 
Local Group dwarf IC 10). 

From these images, we created optical and near-IR 
CMDs of the stars in KKH 98, and estimated the star 
formation and metal enrichment history of the galaxy. 
On average, KKH 98 has experienced a low level (~ 
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5 x 1O~ 4 M yr^ 1 ) star formation rate over its lifetime, 
and has maintained a metallicity of roughly l/30th - 
l/10th solar. In addition, there is evidence for a pro- 
longed period of quiescence (1-3 Gyrs ago) where lit- 
tle star formation occurred. Despite the relatively small 
number of IR bright stars (592 compared to 2539 in the 
optical), the results from the IR data match those from 
the optical for the bulk of cosmic time. 

However, there are some discrepancies in the optical 
and IR derived SFHs for ages < 1 Gyr. In the optical 
CMD, younger populations were primarily constrained 
by the main sequence turn-offs, whereas the IR CMD 
constrained younger populations via the asymptotic gi- 
ant branch. Discrepancies between the two can be ex- 
plained by an over-production of AGB stars in the SFH 
models. We estimate that the models over-predict the 
numbers of AGB stars in KKH 98 by as much as a fac- 
tor of three. The over-production of AGB stars in SFH 
models has bee n noted previously in stu dies of low metal- 
licity galaxies (jGullieuszik et al.l 120081 ) . suggesting that 
AGB lifetimes (and mass loss rates) may be metallicity 
dependent. 

Because the IR bright stars were not crowded, these 
techniques may be useful for measuring stellar popula- 
tions at much larger distances, especially with the James 
Webb Space Telescope and AO on future Thirty-Meter 
class telescopes. However, improved models of AGB evo- 
lution may be required to take full advantage of these 
techniques in the IR. 
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TABLE 1 

Summary of MATCH derived results 



Optical CMD 



IR CMD 



<SFR> M yr~ 
< M/H > 
Reddening (Ay) 
Distance 



.-i 



5.18 x KT 4 

-1.33 
0.52 ± 0.06 
27.06 ±0.09 



5.45 x 10~ 4 

-1.21 
0.41 ± 0.08 
27.10 ±0.07 



TABLE 2 

Stellar density in different evolutionary phases of the CMD 
for both the KKH 98 data and MATCH models. 



Optical CMD Optical Model IR CMD IR Model 



a F814W < 23.2 and F475W - F8UW > 1.5 
b K < 21.2 

c 23.2 < F8UW < 25.5 and 1.5 < F475W - F8UW < 2.5 

d 21.2 < K < 22.5 and 1.0 < F814W - K < 2.5 

c 27.0 < F8UW < 23.0 and F475W - F8UW < 0.5 



AGB 

upper RGB 
MS 



42 a 
851 c 
446 c 



130 
800 
421 



41 b 
207 d 



92 
191 



